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ABSTRACT 



An algorithm for predicting the Maximum Usable Frequency 
between two points on the surface of the earth for High 
Freauenca communications is implemented on a TI-59 hand-held 
programmable calculator. 



1. INTRODUCTION 



An algorithm for predicting the Maximum Usable 
Freauency (MUF) in the High Freouency radio band has been 
developed at the Naval Ocean Systems Center* San Diedo* 
California* by R.B. Rose and J.N. Martin II 1 3 * and has been 
implemented in both FORTRAN IV 3nd BASIC computer languages 
under the name MINIMUF-3.5. The algorithm has been 
incorporated into the CLASSIC PROPHET propagation analysis 
system* and has been demonstrated to provide a 
field-deployable capability for computing HF propagation 
paths with micro-computer resources. 

2. OBJECTIVE 

The essential portions of MINIMUF-3.5 are contained in 
approximately 80 BASIC language statements* reproduced in 
appendix A* which suggests a potential for being encoded in 
hand-held programmable calculators. H.F. Hite of Hughes 
Aircraft Company adapted a restricted version of MINIMUF-3.0 
for the HP-67 calculator that computes MUFs for 1-hop 
F-layer propagation and is thus limited to station 
separations of 4000 kilometers or less. A complete 
implementation of MINIMUF* capable of 2 -hop path prediction* 
was desired. 
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3. PROGRAM DEVELOPMENT 



A TI-59 calculator program was adapted directly from 
the Rose and Martin BASIC program listing with some logic 
changes required to accomodate the memory constraints of the 
calculator. The program listing* summary of logic changes* 
operating instructions* and test case are contained in 
appendices B* C* D* and E* respectively. The program listing 
is annotated with cross-references to statements in the 
BASIC listing for assistance in tracing logic flow. The bulk 
of the development task was to minimize program and data 
storage reoui rements ♦ Several iterations finally produced 
the current version which reauires 800 program steps and 20 
data registers* the full capacity of the normal ly-accessible 
calculator attributes. 



4. PERFORMANCE 

The TI-59 program was optimized for storage* that being 
the overriding constraint. Execution time for single— hop 
predictions is approximately 50 seconds. Were the program 
optimized for speed* given a larger storage capacity* 
execution time could conceivably approach 40 seconds. The 
test case contained in appendix E is the identical test case 
promulgated with the Rose and Martin report* and yields the 
same results on the TI-59. Notice? however* that the MUF 
limit of 32 MHz has been raised to 50 MHz* thereby providing 
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for HF path prediction during periods of high solar flux 



when customary band limitations are exceeded. The test case 
should be executed after initially programming the TI-59 by 
keystroke to ensure correct program entry. A copy of the 
program is available on magnetic cards by sending two blank 
cards to the authors* however* the user is cautioned that 
magnetic cards are not guaranteed to be transportable among 
all TI-59 calculators. 



5. CONCLUSIONS AND RECOMMENDATIONS 



A capability for predicting Maximum Usable Frequencies 
in the HF radio band can be realized with hand-held 
calculators* 3nd therefore can be widely deployed in the 
field. An even larger capability can be realised by 
utilizing custom made modules for the TI-59. Each module 
provides for 5000 program steps in addition to the normal 
960 maximum step capacity of the calculator* and could 
provide for additional portions of the CLASSIC PROPHET 
System to be included in calculator implementations* such as 
the D-region absorption model. Appendix F lists several 
persons to contact in regard to having custom modules 
manufactured for the TI-59. 



APPENDIX A 



BASIC MINIMUF-3 ♦ 5 PROGRAM LISTING 



The BASIC 1 andua^e listing of MINIMUF-3. 5 is reprinted 
with the permission of R.B. Rose* and serves as a guide 
the TI-59 keystroke listing. 



here 

for 
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MlMMLFo.5 FIUX. RAM 



1060 REX 

1018 K7»SIN<L1>5SIN<L2>+C0S<L1>4C0S(L2>*C0S<W2-W1 > 

1029 G1*ACS<X7 MAX -1+1.0E-3 MIN 1-1.0E-8) 

1936 K6-1.S91G1 
1048 K6*K6 MAX 1 
1050 K5*l/K6 
I860 J9*188 

1878 FOR X1«1/<2»X6> TO 1-1/<2*K6> STEP 8.9999-1/K6 
1088 IF K3*l THEN 1100 
1090 K3-0.5 
1100 P*SIN<L2) 

1119 Q*C0S<L2> 

1129 A*<SIH<Li >-P*COS<Gl ) )/<Q*SIN<Gl > ) 

1138 B*G1$X1 

1148 C*P3C0S < B ) +Q*S I N < B ) $ A 

1150 D*<COS<B)-C*P>/(Q$SQR< 1-C$C) > 

1160 D»ACS<D MAX -1+1.0E-5 MIN 1-1.8E-8) 

1170 H6*M2+SGN<SIN<M1-W2) )*D 
1188 IF M0*>8 THEN 1200 
1196 H0-W0+P1 
1208 IF H8<P1 THEN 1220 
1210 H9«H0-P1 

1229 L9«P9-ACS<C MAX -1+1.0E-5 MIN 1-1.0E-8) 

1238 Yl*8.0172*<18+<M0-m38.4+D6> 

1240 Y2*8.489SC0S<Y1> 

1258 X8*3. 822M8+12+8. 13$<SIN< Y1 >+l . 2*SIN<2*Y1 > > 

1269 K8*X8-12*< l+SGN<K8-24) )*SGN<ABS<K8-24) > 

1270 IF COS<L0+Y2) >-0.26 THEN 1330 
1280 K9*0 

1298 G8*0 
1300 M9*2.5*G1*X3 
1310 M9-M9 MIN P9 
1329 M9*SIN<M9) 

1338 M9*1+2.S$H9*SQR<N9> 

1348 GO TO 1599 

1350 X9-<-0.26*SIN<Y2>*SINa.0>>/<COS<Y2>*CQS<L0>+1.0E-3> 

1360 X9«12-ATN<K9/SQR< ABS< 1-K9*K9> > >*7. 639437 

1379 T*K8-K9/2-H2*<l-SGN<X9-K9/2>>*SGN<A8S<K8-K9/2>> 

1388 T4-K8+X9/2-12$<l+SGN<X8+K9/2-24))$SGN<ABS<K8+X9/2-24)) 
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1390 C0 s ABS<COS(L0+Y2> > 

1400 T9*9.7*C3t9.6 

1410 IF T9 >0. 1 THEN 1430 

1420 T9=0. 1 

1430 M9=2.3SG1*K3 

1440 M9=M9 MIN P0 

1450 N9*SIN<H9> 

1468 H9«H-2.3*H9*SQR<M9> 

1470 IF T4<T THEN 1309 

1480 IF <T3-T>*<T4-T3>>0 THEN 1310 

1490 GO TO 1648 

1500 IF <T3-T4>*<T-T5>>0 THEN 1648 

1519 T6-T5+123< 1+SGN<T-T5> >$SGN< ABS<T-T5) > 

1529 G9*PI$<T6-T>/X9 
1538 G8=*PI$T9/K9 
1548 U*<T-T6>/T9 

1550 G8 ss C05<SIN<G9>+G8t<EXP<U)-COS<G9) ) )/< 1+G8SG8) 

1568 G7*C9$<G8*<EXP<-K9/T9>+imEXP<<K9-24>/2>/<H’G8*G8> 
1570 IF G0*>G7 THEN 1390 
1580 G0*G7 

1590 G2*< 1+S9/250>*H9*SQR<64-58*SQR<G8> ) 

1680 G2-G2SC 1-0. 1*EXP< <K9-24>/3> > 

1610 G2*G2$<l-Kl-$GN<Ll)*SGN<L2m0. 1> 

1629 G2*G2$< 1-0. ISC 1+SGN <A8S<SIN<L0 ) )-COS<L0> ) ) ) 

1638 GO TO 1788 

1648 T6*T5+12S< 1+SGN<T4-T5) )SSGN< ABSCT4-T5) ) 

1650 G8-PIST9/K9 
1660 U*<T4-T6)/2 
1670 Ul*-K9/T9 

1688 G0»C8*<G8*<EXP<U1>+1>>*EXP<U>/<1+G8*G8> 

1698 GO TO 1590 

1708 IF G2> J9 THEN 1720 

1718 J9*G2 

1729 NEXT K1 

1730 J9-J9 MAX 2 MIH 32 
1740 RETURN 
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APPENDIX B 



T 1-59 PROGRAM LISTING 



The keystroke listing of the TI-59 program .implementation of 
MINIMUF-3.5 follows. Segments of TI-59 code are 
cross-referenced to the BASIC program listing. 
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1080,90 



000 


6 


LBL 


050 


•2‘ DL 


x : t 


0 0 l 


I 


H 


051 


6 1 


i 


002 


- 

*r :C 


STD 


052 


6 7 


TQ 


0 0 3 


i y 


18 


053 


0 0 


00 


0 0 4 


Q 9 


PRT 


054 


57 


■J ( 


005 


25 


CLR 


055 


93 




00c 


8 0 


deg f 015720 


056 


05 


j 


G 0 7 


a 3 


RCL 


057 


65 




008 


15 


15 


058 


4 3 


RCL 


009 


— * cr 

. .j 


- 


059 


0 6 


J 6 


0 1 0 


4 3 


RCL 


060 


65 




0 1 1 


13 


13 


06 1 


05 


j 


0 1 2 


95 


— 


U 6 2 


95 


— 


0 1 3 


39 


CDS 


0 6 3 




X i T 


014 


65 


X 


064 


3 9 


it 


0 1 5 


'—r .j: 


RCL 


065 


. 1 1 


S3R 


0 1 6 


i 2 


12 


0 6 6 


07 


0 7 


0 1 7 


9 


CDS 


067 


94 


94 


0 1 8 


65 




0 6 8 


LI* C 
.J. 


— 


0 1 9 


4 3 


RCL 


069 


02 


■ j 


020 


i 4 


14 


07 0 


9 5 


— 


02 1 




CDS 


07 1 




SIN 


022 


85 


+ 


072 


-5 


Y' : 


n 2 3 


. - T : 


RCL 


073 


HI 


1 


024 




12 


074 


M ;T; 




025 


•- O 


S I H 


075 


U 5 


r; 


0 2 8 


i 5 


X 


078 


65 




027 




RCL 


077 


02 


j 


028 


i 4 


14 


078 


9 3 




029 


■ 5 C : 


S I N 


0 7 9 


05 


~: 


030 


95 


- 


080 


35 


f 


03 1 


70 


ROD 


08 1 


0 1 


1 



032 


• j 


1 HV 


032 


39 


CDS 


0 3 4 


4 2 


STD 


035 


. i * 6 




036 


if. 5 


x 1030,40,50 


037 


01 


1 


0 3 8 


93 




0 3 9 


05 




040 


09 


9 


04 i 


95 


= 


042 


u. 


j- 

: • i 1 


04 3 


01 


1 


04 4 


"1 


S 3 R 


045 


ii7 


07 


04 6 




95 


047 


85 


1 XX 


048 


42 


era 


0 4 8 


07 


07 



G 82 

083 

084 

085 
083 

087 

088 

089 

090 
09! 
0 9 2 

093 

094 
0 95 
09 6 
0 97 
098 
0 9.9 



95 

42 STD 
08 08 
01 i 
0 0 0 
00 J 

■-2 STO 

1 11 
c . _ _ 



“3 RCL 
or 07 
09 INI 
95 - 

42 STD 

U 9 o 9 
i 1 i l 



1300,10,20,30 



1060 



1070 



1070 - cont'd 



10 



100 


43 RCL 






150 


5 




1 0 1 


ij f U f 






1 5 1 


--- . -J 


R C L 


102 


95 






152 


i i 0 


00 


103 


R *< 






■f nr •-» 


8 


S I N 


1 04 


y 3 ( 






154 


*5 




105 


43 POL 






155 




i< 0 L 


1 00 


0 9 0 9 






156 


01 


01 


1 0 ? 


cr 

•J 






157 


■: n~ 
J 




1 0 £ 


0 1 1 






158 


•-! 


RCL 


109 


54 > 






153 


: 2 




1 10 


If; 5 + 






1 60 


;-s 5 




1 1 i 


4 3 RCL 






1 0 1 


•- 3 


L 


112 


0 7 0 7 






162 


i i 0 


0 ij 


1 1 3 


55 4 






163 


; ; 9 


C 2 S 


114 


02 2 






164 


^5 


■“ 


1 15 


0 5 






1 65 




S T l 


1 1 6 


ijO X 


1130 




166 




j 3 


2 - 5* 


‘43 RCL 






167 




I My 


1 i S 


i j 6 0 8 






lc-S 


■ s 


CDS 


119 


3 5 






169 


: '■ 0 




120 


42 8 TO 






170 


" • 5 


t- 


121 


00 GO 




- 


1 T’ 1 


: 9 


n 


i ill i 


43 RCL 


1100 ,10,20 




172 


’ D 


- 


1 23 


4 1 4 






173 


: 2 


2 


124 


!-. 0 DEG 






174 




-■ 


125 


39 COS 






175 


30 


STD 


126 


42 8 TO 






176 


f! 


n 


12 "' 


0 1 0 1 






177 


- 3 


RCL 


128 


35 1/X 






173 


0 0 


0 0 


129 


65 x 






179 


' ; Q 


c : J 3 


1 30 


53 <: 






1 80 


’ “ 


- 


1 3 1 


43 RCL 






1 8 1 




RCL 


132 


12 12 






182 


n 3 


u 3 


133 


38 SIN 






183 


1 .5 




134 


’’5 






1 34 


-• *7: 


RCL 


135 


43 RCL 






135 


ijV 


2 


136 


4 14 






1 8 6 


: : * 


- 


1 37 


33 SIN 






1 8 


c 

... .j 


- 


138 


‘-2 STQ 






188 


•Vfl 


R .'L 


139 


i ! 3 y 3 






1 82 




•J 1 


1-10 


6.5 






1 4 u 


-5 


- 


3 41 


4 3 RCL 






191 


i 1 




142 


n 8 06 






192 


M] 


j 


1 43 


’0 R 4 D 






1 9 3 


• cr 

_.i 


- 


1 44 


39 CO 8 






194 


• -31 


R L. 


i =• 


95 - 






145 


■" 1 


J '.~i 


I 4 r. 


u- cr 

. i ...1 T 






1 96 


3 


X -* 


•147 


43 RCL 






i9 T 


* • j 


j 


1 48 


06 ’ 06 






198 


j. 


r m 


1 4 9 


SIN 






199 


<5 


= 



1140 



1220 



1150 



11 



200 


2 2 


I H V 


1160,70,80,90,1200,10 


250 


■7 


17 


201 


•0 9 


cos 




251 


. : cr 

.J 




202 


6 5 






252 


!j 3 


■ 2 - 


203 


C 7 •7: 


(* 




253 


uo 


J 


2Q- 


-8 


ROL 




254 


4 3 


, 


205 


j. - 2 = 


13 




255 


04 


-• 


206 


' ’5 


- 




2 0 4: 


-5 




207 


-3 


ROL 




257 


53 


< 


208 


. 5 


15 




nr 

21 -J C : 


43 


ROL 


209 


54 


} 




259 


6 


16 


2 1 0 


6 0 


DEG 




230 


•5 


- 


211 




S I N 




261 


0 1 


1 


2 1 2 


6 9 


OP 




2 6 2 


95 


nr 


“i ■>. 


■ 0 


10 




2 6 3 


'-r nl 


3 TU 


21 d 


8 5 


4- 




264 


nn 


00 


2 ] 5 


-T .S 


ROL 




265 


70 


R H jj 


2 1 6 


■ nr 


15 




2 6 6 


■3 9 


C 2 S 


217 


6 5 






267 


05 




2 10 


8 9 


n 




263 


49 




2 i 9 


r “ i 


- 




2 6 9 


1 j 4 


•4 


220 


ill 


i 




270 


0 0 




221 


0 8 


3 




271 


n 9 


9 


2 2 2 


00 






■l i. f 2 


95 


::: 


”■ I s 

J_ i— 


r D 


= 




2 ?' 0 = 




370 


224 


''j *7* 


X 1 1 




274 


0 1 


j 1 




8 9 






2 •’ 5 




X : T 


4 


b5 






276 


“ 

“s . J. 




-7 


: 7 : 








■: 








.2. 




2 = f 


• J J: 


2- 


_ sl_ 


!T* 


4 




2 s' 0 


9 3 


;1 


229 


O 2 


K ; T 




279 


; j 3 


7 ; 


230 


9 5 


- 




2 8 0 


02 




231 




IN V 




261 


. 

.J 


4 


'I s *7* *3 


•' 7 


OE 




2 0 2 


01 


1 


2 0 3 


02 


02 




2 0 1= 


02 


2 


234 


•j C : 


33 




234 


05 


+ 


nr 

p . j 


— • rr 


- 




235 






2 3 6 


* 2 


< i l 




2 8 6 


i'U 


1 


iL • "j :* 


9 5 


:z: 




7 


i 1 3 


“ 1 




.. 


•' > T 




7 7 ; L : 


nr 




“ *•"* •-* 


•-* l' 


. i • 




2 ; i ; 0 


•: : s 




2 39 


93 


s 


1230 


Q 

1 — -• * 


r .3 


< 


240 


00 


;J 




2 9 0 


-3 


ROL 


241 


0 1 


1 




291 


uo 


0 0 


242 


07 


s' 




2 9 2 


■ •; P 


Sli-I 


243 


02 


:• 




2 9 3 


5 




2 4 4 


f:,5 






294 


Ui 


I 


245 


C* -1 

• 1 0 


i’ 




2 9 5 


! ? 3 




24 6 


0 : 


1 




296 


i i 2 




247 


00 


0 




297 


6 5 


X 


24 8 


* •. c* 
1 1 - 1 


f 




2 9 3 


l; .-i 

. > j 


f 


2 4 9 


43 F 


:cl 




2 9 9 


43 


ROL 



1240 



1250 



12 



! 0 0 


o n 


Ci 0 


350 


4 


STO 


:01 


65 




351 


U 3 


03 


1 1 J V 


Ij2 


— 


352 


•0 1 ~ 


r x 


:03 


•54 




x 5 J 


65 




: 0 4 


X 3 


S I N 


354 


05 


rr 


05 


•45 




355 


0 3 




:06 


-■* ir 
J 


- 1260 


356 


«. i -r 


-h 


107 


3 2 


x ; i 


357 


i j 6 


6 


: u c 


ijv 




£T 

•0 : • J c* 


95 


- 


;Q9 


04 


4 


359 


34 


fx 


' i 0 


8 ?5 


= 


360 


65 




HI 


3 2 


x ; t 


361 


4 >! 


r :l 


! 12 


CT 

.j. 


- 


3 6 2 


03 


0 8 


•13 


01 


i 


363 


6 5 


x 


H4 


02 




364 


IT ~ 


: 


H 5 


05 




365 


0 i 


1 


16 


53 


if 


X b X. 


85 


+ 


1 7 


01 


i 


3 6 7 


4 3 


RS L 


1 1 b 


35 


4 


3 6 8 


A 9 


15 


1 9 




X i I 


X b 


r- r 
• _• • J 


4. 


20 


0 9 


08 


370 


02 




21 


‘ 0 


1 0 


371 


05 


- 




4 2 


STO 




00 


6 




uo 


JIJ 


• • -Z : 


95 




4 S T 


54 


) 


374 


65 




sr 


6 5 




375 


L 

• 0 : 


•' 


\- s ■ 


4 3 


ROL 


3 R 6 


0 I 


i 


27 


0 0 


0 0 


•Z : : ; 


* ? 5 


_ 




50 


I x I 




*7 3 




.29 


95 


-- 


3 ■*’ 3 


0 1 


- 


30 
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APPENDIX C 



PROGRAM LOGIC MODIFICATIONS 



The following program logic changes were made to the BASIC 
version of MINIMUF-3.5 to minimize program storage 
reoui remerrts t 

(1) The computation of M9 at lines 1300-1330 and 1430-1460 
is independent of the main loop and W3S moved to the 
beginning of the program* immediately following the 
calculation of G1 and K5. 

(2) The loop index computation at line 1070 was separated 
into a computation of a new variable* KHGP* which is either 
.1 or 2* and a follow-on computation of K1 from KH0P as the 
first item within the loop. 

(3) The computation of GS at lines 1530 and 1650 is common 
to both legs of the preceedinS test for (T5-T4) (T-T5) * and 
has been moved ahead of that test. 

(4) The test at line 1470 serves to reverse the sense of the 
following tests at lines 1480 and 1500. This logic has been 
combined into a single test on ( T5-T4 ) ( T-T5 ) with a 
weighting factor of 1 or -1 to reverse the sense of that 
test ♦ 
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(5) The computation at line 1360 was simplified to the 



arcsine by way of trigonometric identity* 

(6) The MUF limit established at line 1730 was changed to 50 
vice 32 to accomodate hidh solar flux densities* 

(7) Throughout the prod ram ? explicit calculation and storage 
of variables that are used only once in following statements 
has been eliminated to conserve on data register 
reaui rements ♦ 
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APPENDIX D 



PROGRAM OPERATING INSTRUCTIONS 



The following instructions must be followed to operate the 
T I -59 version of MINIMUF-3,5! 

(1) Repartition the calculator for 800 program steps and 20 
data registers by pressing 2/0P/17. 

(2) Load the 4 memory banks from program cards (2) or by 
entering program keystrokes with the calculator. in the LRN 
mode . 

(3) Enter input d3ta as follows^ 

(a) Transmitter North Latitude in decimal degrees in 
R12 (range -90 degrees to 90 degrees) 

(b) Transmitter West Lonstitude in decimal degrees in 
R13 (range 0 to 360 degrees) 

(c) Receiver North latitude in decimal degrees in R14 
(range -90 to 90 degrees) 

(d) Receiver West Lonstitude in decimal degrees in R15 
( range 0 to 360 degrees) 
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(e) Month in R16 ( range 1 to 12) 



(f) Day in R17 (range 1 to 31) 

(g) Sunspot Number in R19 (range - positive number) 

(4) Enter Time in X-register (range 0 to 24 decimal hours) 

(5) Press A 

(6) If a printer is attached? the input time is echoed? 
followed by the answer (MUF) after 50 seconds for sinS.le-hop 
and 100 seconds for two-hop predictions. 

(7) The answer is displayed in the X-redister when the 
program halts. 

(8) All input Quantities remain undisturbed in R'12 thru R19 
(time is stored by the program in R13). Another MUF for a 
different time can be computed directly by repeating steps 
(4) and (5). 
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APPENDIX E 



TEST CASE 



The test case that follows was provided in the NOSC report 
on MINIMUF-3.5. The actual printer listind of input and 
program output is included. 






CONTENTS 
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INPUT 
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APPENDIX F 



TI-59 CUSTOM MODULES 



The following points of contact are provided for the design 
and manufacture of TI-59 custom modules. Both Texas 
Instruments and Horizons Technology (contracting with TI) 
provide software? emulator? and consulting support for the 
design and Production of custom modules. Costs Quoted at 
this time are approximately $12r000 for a minimum order of 
250 modules. This includes emulator support and manufacturer 
but does not include software consulting. Cost reduction 
through Quantity is available. 



Fred Wilke 

TI-59 Custom Module Division 
Texas Instruments 
Lubbock r TX 
(806) 741-3240 

Robert Kruger 

Horizons Technology? Inc. 

7830 Clairemont Mesa Boulevard 
San Diego ? CA 92111 
(714) 292-8331 
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